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Abstract 

“'“n waves in deep water interact with a current, the direction of propagation 

ocean “ ™ such “ hei « hl “ d length are affected. Swell in the open 

Mem an undergo significant refraction as it passes through major current systems like 
the Gidf Stream or Antarctic Circumpolar Current. Remote sensing techniques such as 
synthetic aperture radars (SAR) have the potential to detect wave systems Z a ~d“ 
geograplncal area. Combining a model for wave refraction in the presence of currents 

sdL PIOblein ° f ,he wave data can be 

° , d ‘° determine the direction and magnitude of the intervening currents In this 

is erUi Tl" o SWdl meaSUIed br SAE <® a stellite pass over the Gulf Stream 
Tbe predicted by a numerical model under conditions of varying 

current profiles mid velocities is compared to SAR generated wave spectra. By matching 

dtla^T! P ^ IT “ ‘ he bSt correlation ^ve refraction to the SAR 
datiq the tomographic problem of measuring the Gulf Stream current is solved 

is mS!l COTI . r 1 b T tW ' en ‘ he m ° dd “ d SAR data “ ob'sined when a current 
2 Th A a ° P h ’ ,dOClt7 Pr ° Me With a dkection of and a current speed of 
2 m/s. The direction agrees with that visually observed from the SAR images, and the 

ZoTZcpsTat "* d “ e t0 ,ie Coas ‘ fruard estimates for the Gulf Stream at the 
ime of the SEASAT pass. The current profiles used did not take into account a possible 
widening of the Gulf Stream at the position of the satellite overpass. There is a great 
deal of scatter m the SAR data, both before and in the Gulf Stream, so it is diffiodt to 
correlate every point with specific current behavior, but the increase in wave length and 
change m wave angle m the center of the Gulf Stream seem to indicate that them may 
be a non-uniform feature such as the formation of an eddy or other lateral variability 
near tne current s edge. J 
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Chapter 1 



Introduction 



1.1 Motivation and Problem Statement 

When ocean waves in deep water interact with a current, the direction of propagation 
and characteristics of the waves such as height and length are affected. The influence of 
currents on wave propagation is more pronounced in deep water, because as the waves 
shoal, depth refraction dominates the changing wave characteristics. Swell in the open 
ocean can undergo significant refraction as it passes through major current systems like 
the Gulf Stream or Antarctic Circumpolar Current. Remote sensing techniques such as 
synthetic aperture radars (SAR) have the potential to detect wave systems over a wide 
geographical area. Combining a model for wave refraction in the presence of currents, 
with SAR. measurements, the inverse problem of using the measured wave data cam be 
solved to determine the direction and magnitude of the intervening currents. In this 
study the behavior of swell measured by SAR on a satellite pass over the Gulf Stream 
is ex amin ed- The refraction predicted by a numerical model under conditions of varying 
current profiles and velocities is compared to SAR generated wave spectra. Complex 
correlations are calculated between the vectors formed by SAR measured wave length 
and direction and the wave length and direction vectors predicted by the numerical 
model. By chosing the model inputs which result in the best correlation, the tomographic 
problem of measuring the Gulf Stream current is solved. 
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In a pioneering work, Unna (1942) showed how waves axe shortened or lengthened 
by opposing or following currents, respectively. Johnson (1947) derived a form of Snell’s 
law which describes the refraction of wave crests when interacting with a current at an 
oblique angle to the wave propagation direction. Arthur (1950) introduced the concept of 
wave energy transmission along rays which are not necessarily perpendicular to the wave 
crests. Longue t- Higgins and Stewart (1960) demonstrated the importance of nonlinear 
interactions between waves and currents and introduced the concept of radiation stress to 
describe the variation of wave energy in the presence of a moving medium. Kenyon (1971) 
derived an analytical solution for the refraction of orthogonals and rays in the presence 
of a linear current shear. Abemethy and Gilbert (1975) developed a numerical model to 
predict the refraction of surface gravity waves over an ocean bottom of varying depth. 
Treloar (1985) modified this model to include the linear effects of current interactions. 

A depth refraction model developed by Brampton (1977) is modified to solve the 
kinematic problem of wave refraction in currents. This current refraction model is val- 
idated by comparing it to analytical solutions derived by Johnson and Kenyon. It is 
then applied to an ocean situation, and the wave kinematics observed by a synthetic 
aperture radar flown on SEASAT during a pass over the Gulf Stream are compared to 
the kinematic behavior predicted by the model. 



1.2 Outline of Present Study 

In section 2 the theory of linear wave interactions with currents is reviewed. The 
governin g equations of wave refraction in the presence of currents are shown, and the 
derivations of analytical solutions to the refraction of wave fronts passing over a disconti- 
nuity, and waves passing into a linear shear current are reviewed. Then the assumptions 
and structure of the numerical model used in this study are described, and a number 
of examples are shown comparing analytical and numerical solutions to the problem of 
wave refraction. 

In section 3 the principles of imaging real waves with synthetic aperture radar are re- 
viewed, and the characteristics of the instr um ent as deployed on SEASAT are described. 
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The SEASAT data from pass 1339 used in this study are discussed. 

In section 4 the application of SAR data and the numerical model for predicting 
wave refraction in the presence of a current is applied to the Gulf Stream. First the 
environmental and meteorological conditions which existed at the time of pass 1339 are 
presented. Next the analysis of SAR images to produce wave energy spectra is described. 
Using these spectra, a wave system’s propagation is traced across the Gulf Stream. This 
is compared to the predicted path of wave rays and orthogonals over the Gulf Stream. 
A number of different possible current profiles for the Gulf Stream are examined. By 
correlating the predicted wave ray paths and wave lengths with SEASAT measured 
spectra, a best fit current profile is proposed. 
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Chapter 2 



Theoretical Background 



2.1 Theory of Wave- Current Interactions 

2.1.1 Governing Equations 

The study of the kinematics of wave-current refraction really involves two problems: 
(1) the direction of wave propagation, and (2) the changes in wave length. The wave or- 
thogonals (perpendicular to wave crests) follow paths tangent to the wave number vector 
at a phase speed modified by current interactions. This is the direction of propagation 
seen by an external observer who takes a ’snapshot’ of the wave system. The wave ray 
gives the direction of energy propagation traveling at the group velocity and is given by 
the vector sum of the intrinsic group velocity and the current velocity vectors. 

The assumptions made in this study in applying the governing equations for wave- 
current interactions are as follows: 

1. The effects of viscosity and surface tension are neglected. 

2. The fluid is assumed to be incompressible and its motion nearly irrotational. 

3. The currents are assumed to be strong and vertically uniform. 

4. The spatial and temporal variations of the horizontal currents are assumed to be 
slow relative to the wave scales; the fluid is deep relative to the wave length (short- 
wave approximation). 
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Figure 2-1: Wave ray and orthogonal relationships 

5. Reflection and diffraction effects are excluded. 

6. The wave amplitude is assumed to be small. 

7. Effects by the surface curvature of the earth axe small enough to be neglected. 

For a steady, inhomogeneous medium, the wave frequency is constant and hence the 
absolute phase speed (c tt ) can be found from the relative or intrinsic frequency (o> r ), 
wave number vector k, and current vector U, i.e., 



u r + k • U 


(2.1) 


Cr + U sin a 


(2.2) 



where a is the angle between the wave orthogonal (or the wave number vector k) and 
the current normal, as defined in Figure 2-1. 

The path of wave orthogonals is determined from integration of the equations: 



dx 

ds 


= cos 6 


(2.3) 


dy 

ds 


= sin0 


(2.4) 


dd 

ds 


sinfl&f-cosO^ 

C a 


(2.5) 
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where Ca — ^ is the absolute phase speed, as seen by a stationary observer, 6 is the 
angle between the orthogonal and the x-axis, and an elemental distance along the arc 
is given by ds = c a dt. The relationship between wave crest and current is illustrated in 
Figure 2-1. 

The wave ray path is determined from integration of the ray equations (Kenyon, 
1971): 



dx du 

~dt = dk 
dk _ du 
dt dx 



( 2 . 6 ) 

(2.7) 



Here is the group velocity; describes the change in wave number along rays. 
Another form of the ray equations (Treloar, 1985) is: 



dx 

dr 

dy 

dr 

d[ i 

dr 



cosn 



sin h 



cos 2 h 



d t an fi 
dr 



( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 



where r is now the distance along the ray, and the angle /z is defined by: 

Cgr sin 9 4- v 
tan h = — — 



( 2 . 11 ) 



T COS 0 d" XL 

where is the relative group velocity, 9 is the direction of the wave orthogonal, and u 
and v are the x and y components of current U. The relationship between the current 
and wave ray is also shown in Figure 2-1. 



2.1.2 Analytical Solutions 

Two analytical solutions are described here to demonstrate the applicability of these 
equations. These solutions are used to validate the numerical method outlined in 2.1.3. 



Johnson’s Derivation of Snell’s Law 

Johnson (1947) derived the path of wave orthogonals entering a current from still wa- 
ter. This situation is illustrated in Figure 2-2. The directions of wave crest propagation 
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Figure 2-2: Wave orthogonal entering a constant current from still water 



sin p 



and currents are related by the following expression: 

C Q 

sin a 

The initial or still water phase velocity is Cq, and the phase velocity relative to the 
current is c. The waves must be continuous across the discontinuity, so the wave length 
in still water L a is related to the wave length in the current L: 

L _ L 0 
sin/3 sin a 

Equation (2.12) can be expressed as: 

csina 



( 2 . 12 ) 



(2.13) 



sin/3 = 



(2.14) 



Co - U sin a 

Because in deep water phase speed is related to wave length by c 2 = |^, where g is the 
gravitational constant, (2.13) can be rewritten as: 
c 2 c 2 

4-j = 4 s - (2.15) 

sm p sin a 

A form of Snell’s law can be derived from (2.14) and (2.15). Upon simplifying: 

sin a 



sin/3 = 



(2.16) 



' (1 — ^ sin a) 3 

The final result (2.16) describes the refraction of wave crests from still water across a 
current discontinuity. 
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Figure 2-3: Wave crest entering a current which has a constant shear 

Kenyon’s Solution for the Case of a Linear Shear Current 

Kenyon (1971) used the ray equations to derive an analytical solution to both ray 
end orthogonal paths for the case of a current with a linear shear. The ray equations 
hold under the geometrical optics approximation that wave amplitude and frequency 
vary slowly over distances of the order of the wave length. 

An approximate equation for the radius of curvature of the rays is: 

R = ^ (2.17) 

where £ = - fjj), which is the component of vorticity in the positive z direction. 

The current velocity is assumed to be in the x-direction and only a function of y; it 
has a constant shear a: 

The initial conditions are that the rays pass through * = 0, y = 0 at time t = 0 with 
initial angle <f> 0 = cot -1 jf 2 . This situation is shown in Figure 2-3. 

The integration of the ray equations (2.6-2. 7) using the above initial conditions results 



ay if y > 0 
0 if y < 0 



(2.18) 
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MO = t 

y(f') = 2 i 2 <,(siii<i 0 ) - 2 [(1 + t 2 )* — (sino 0 ) - j] 



(2.19) 

( 2 . 20 ) 



in: 



z(t) = + 

cot<? 0 (sinp 0 )-i +2->[f , (0 > iT)-f’(© o> ir)]} (2.21) 

where t = cot 6 0 — st , F(Q,K) is the elliptic integral of the first kind, cos 0 = (l-f-t <2 ) - « , 
cos © 0 = (sinp 0 )i, and K = 2"*. In (2.20) and (2.21) R* = This is the exact 

initial radius of curvature of the rays, where im q is the initial wave frequency. 

Kenyon also showed that Johnson’s form of Snell’s law can be put in the form x = x{y) 
for the same shear current by using tand = 4^. This results in a description of the wave 
orthogonals: 

x (y) = - f(©..J0] (2-22) 

where cos 0 = (sin do)* (1 + y ^f a ) -1 - In Kenyon (1971), 0 is mis defined as cos 2 0 = 
sm^(l + 2ic£*)-i. 

2.1.3 Numerical Solutions 

Over the past four decades numerous graphical and numerical methods have been 
developed for refraction calculations (see for example, Chao (1972) and Earle and Madsen 
(1987)). A somewhat different algorithm has been first described by Pararas- Carayannis 
et aL (1968) and later by Abemethy and Gilbert (1975). Both algorithms are based on 
triangular grid systems which assume that within each triangle the velocity field is well 
approximated by a plane. For this condition, the ray path can be solved analytically 
and in particular, traces out the arc of a circle. There are three major advantages in 
assuming V c constant within a triangular element over the conventional approach with 
a rectangular grid. One, an exact solution of (2.5) is possible for triangles but not for 
rectangular elements. Two, iterative predictor-corrector techniques axe usually required 
to advance the ray in short steps across a rectangular grid element. Three, the ray 
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// 

t 

CENTER OF ORCIJE 

Figure 2-4: Ray path along a constant circular arc through a triangle 

path calculations in triangles are much faster because there is no need to specify a time 
step nor a desired level of accuracy for the iterative procedure. Abemethy and Gilbert 
(1975) used equilateral triangles in their algorithm. An example of the application of this 
algorithm can be found in Graber et al. (1989), where the variability of wave kinematics 
and dynamics due to depth refraction in the North Sea is examined. Brampton (1977) 
modified Abernethy and Gilberts’ algorithm by using right triangles. These elements 
allow variable grid spacing in the two dimensions and simplify the angular relationships 
within the triangular grid. A typical triangular element can be seen in Figure 2-4. 

For this study the algorithm used by Brampton is further modified to include the 
effects of currents in deep water. The current field is represented by discrete current 
vectors at the vertices of triangular elements as seen in Figure 2-4. The orthogonal 
equations (2.3), (2.4), and (2.5) are applied to the elements. Using (2.16) the absolute 
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wave velocity at a point in a current can be calculated from: 

C * = (1-£os 7 )* 

where 7 is the angle between the wave orthogonal and the current shown in Figure 2-1. 
By inserting a current field with the values specified at each comer of the triangular 
element, a new phase speed plane can be calculated which includes the current effects. 
The circular paths traced out correspond to the wave orthogonals. 

The calculation of ray trajectories is more complex, since the orthogonal path across 
each triangle must be computed first. The group velocity at each comer is calculated 
by taking the phase speed resulting from the wave-current interactions. The absolute 
group velocity is approximated by adding the current components normal to and along 
the group velocity vector at each point. In deep water c g = \c. The entering angle of 
the ray, fi, must also be calculated for each element. 






(Cgr + U COS 7) 2 -f (17 sin 7 ) 2 ] 2 
+ tan-( ^2 ) 



(2.24) 

(2.25) 



Each triangle must be analyzed twice, first to find 7 from the orthogonal path and then 
to trace the wave ray. 



2.2 Examples of Solutions 

In order to validate the numerical solution technique described in the previous section, 
a number of scenarios have been posed in which the numerical solutions can be compared 
to analytical results. 

The simple case first addressed by Johnson (1947) involves the passage of wave sys- 
tems from still water across a discontinuity into a steady, -uniform current. In the figures 
in this section, the path taken by wave orthogonals are marked by thin lines and the 
wave rays are shown by thick lines. Table 2.1 compares results from the analytical and 
numerical solutions. 
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(b) 10 sec. wave entering a 2 m/s opposing current 




Figure 2-5: Current refraction examples of 10 second waves entering 2 m/s currents, (a) 
Following current, a = 30°. (b) Opposing current, a = 45°. Note that a heavy line 
represents the ray path and a thin line the orthogonal path. 
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Figure 


Current 
(2 m/s) 


Q 


Analytical 

Orthogonal 


Numerical 

Orthogonal 


Analytical 

Ray 


Numerical 

Ray 


2-5(a) 


Following 


30.° 


34.8° 


35.° 


44.6° 


44.9° 


2-5(b) 


Opposing 


315.° 


323.5° 


323.5.° 


338.4° 


337.8° 


2-6 


Following 


55.° 


Reflected 


Reflected 


Reflected 


Reflected 



Table 2.1: Comparison of analytical and numerical results for waves entering a constant 
current. 

A 10 second wave which enters a 2 m/s following current at an angle greater than 
53° from the current normal will be reflected. This angle is known as the critical angle. 
This can also be seen by the solution generated by numerical method in Figure 2-6. 

A more complicated situation involves a wave system which enters a current with a 
linear shear. Figures 2-7 through 2-9 compare the outputs generated by Kenyon’s ana- 
lytical solution with the numerical solution. It can be seen that the analytical solutions 
are almost identical to the numerical solutions to each problem, validating the results 
generated by the numerical model. 

The numerical solution allows us to examine situations which would involve very 
complicated analytical solutions. For instance Figures 2-10 and 2-ll(a) show the situa- 
tion of a ring shaped current loop with a cos 2 <f> profile. It can be seen that waves may be 
trapped, reflected, or simply propagate through such a current system depending upon 
the angle of incidence. 

Another interesting situation is the refraction of a wave field propagating through 
a large scale eddy such as a warm core ring often found near the Gulf Stream. The 
mod eling of such a warm core ring has been described in Mapp et al. (1985). The flow 
field used to simulate this large eddy was a steady circular jet. The radial profile of 
tangential velocity is a power function joined to a Gaussian. This example can be seen 
in Figure 2-ll(b). Orthogonals which penetrate into the center of the warm core ring 
were refracted, but not greatly. Orthogonals which grazed the left side of the ring were 
reflected; they essentially entered a following current at greater than the critical angle. 
One orthogonal which entered the right side of the ring became trapped in an opposing 
current and was refracted almost 90°. 
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10 sec. wave entering a 2 m/s following current 




Figure 2-6: Example of total reflection for a 10 second wave entering a 2 m/s following 
current, a = 55°. Note that a heavy line represents the ray path and a thin line the 
orthogonal path. 
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1 1 sec. wave propagating through a following current 




1 1 sec wave propagating through a following current 




Figure 2-7: Refraction of 11 second wave entering a following shear current which varies 
from 0 to 2 m/s, a = 60°. (a) Solution generated using Kenyon’s analytic result, (b) 
Solution generated with numerical model. Note that a heavy line represents the ray path 
and a thin line the orthogonal path. 
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(a) 1 1 sec. wave propagating through an opposing current 




1 1 sec. wave propagating through an opposing current 




Figure 2-8: Refraction of 11 second wave entering an opposing shear current which varies 
from 0 to 2 m/s, a = 30°. (a) Solution generated using Kenyon’s analytic result, (b) 
Solution generated with numerical model. Note that a heavy line represents the ray path 
and a thin line the orthogonad path. 
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(a) 1 1 scc - wave propagating through an opposing current 





Figure 2-9: Refraction of 11 second wave entering an opposing shear current which varies 
from 0 to 2 m/s, a = 60°. (a) Solution generated using Kenyon’s analytic result, (b) 
Solution generated with numerical model. Note that a heavy line represents the ray path 
and a thin line the orthogonal path. 
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20 second wave trapped in an opposing current 




10 second wave reflected from a following current 




Figure 2-10: Examples of a cos 2 <f> current profile: (a) Wave trapping (20 second wave), 
(b) Wave reflection (10 second wave). Note that a heavy line represents the ray path 
and a thin line the orthogonal path. 
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5 second wave propagating through an opposing current 




(b) xlO 5 11.3 sec. waves proceeding into circular eddy. 




Figure 2-11: Examples of: (a) 5 second wave propagating through a cos 2 <f> current 
profile, (b) 11.3 second wave propagating through a warm core ring. Note that a heavy 
line represents the ray path and a thin line the orthogonal path. 
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Chapter 3 



SAR Imaging of Ocean Waves 

3.1 SEASAT 

One of the instruments deployed on the oceanographic research satellite SEASAT was 
a synthetic aperture radar. This microwave sensor artificially synthesizes an aperture or 
antenna which is hundreds of meters long in space. Therefore, a fine spatial resolution of 
about 20 m can be obtained even from altitudes of 800 km. The basic imaging response 
of a SAR to the ocean surface is explained by Bragg scattering from short gravity waves. 
In the case of radars in the frequency range of SEASAT SAR (1.27 GHz, or the L- 
Band), the Bragg waves are of the order of 20 cm. The Bragg waves, which are formed 
by local winds, are superimposed on underlying long (greater than fifty meter) waves 
which modulate the Bragg waves. Thus, the Bragg waves allow the SAR imaging of the 
longer waves. This implies that one requirement for SAR imaging of ocean swell is a 
local wind field sufficiently strong to create these waves of about 20 cm length. 

One potential use oceanographers see for SAR is in the collection of directional ocean 
wave spectra. Kinsman (1965) and Phillips (1980) have shown that the directional wave 
energy spectrum is probably the most complete single descriptor of the ocean surface. 
When the only source of these spectra is from local measurements taken from bouys or 
ships, their use to oceanographers who seek to understand huge scale ocean processes 
is limited. Remote sensing using SAR promises to provide a valuable tool to better 
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understand the ocean environment. 

The resolution obtained by SAR consists of range resolution (ground resolution per- 
pendicular to the satellite trajectory) and azimuth resolution (ground resolution in the 
direction of the satellite trajectory). 

The range resolution, 6y , is obtained by the pulse length, 6r, projected onto the 
surface. 



Sy = 



(3.1) 



6r _ cr 
2 sine 2 sine 

where c is the velocity of light, r is the pulse duration, and e is the incidence angle. 

The azimuthal resolution 6x is linked to the accuracy in determining the Doppler 
shift of the backscattered signed. A series of pulses are transmitted by the radar; these 



are sampled along the radar trajectory, the synthesized aperture. For a stationary target, 
the azimuthal resolution only depends on the size of the real ant enn a. 



(3.2) 



where D is the real antenna length. 

The imaging of waves by SAR is influenced by three effects. These are: (1) tilt 
modulation, (2) hydrodynamic interactions, and (3) motion effects (velocity b unching 
and azimuthal image smearing). These are discussed in Monaldo and Lyzenga (1986) 
and in Graber and Kelly (1988). Tilt modulation is a geometric effect which results 
from the different incidence angles of Bragg scattering waves seen by the radar as they 
propagate over longer waves. The hydrodynamic effect is characterized by modulations 
in power caused by changes in the spectral density at the Bragg wave number across the 
profile of the long waves. Hydrodynamic modulations result from long wave-short wave 
interactions which increase the ripple amplitude near swell crests and decrease the ripple 
height near swell troughs. They also come from the distortion of the airflow above the 
sea surface by swell waves. These effects are environmentally dependent and generally 
not well understood. Velocity bunching causes image degradation due to the movement 
of the sea surface resulting from the orbital velocities of long waves. The rising face of 
the wave is shifted on the image in the direction of flight, while the falling face is shifted 
opposite the flight direction. For waves traveling in range, this has no imaging effect. 
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However, waves propagating in the azimuth direction are affected significantly. S imil ar 
effects caused by waves shorter than the azimuthal resolution dimension cause Doppler 
smearing. 

The combination of these two effects creates a limitation faced by all SARs viewing 
the ocean: a degradation of the wave energy spectrum in the azimuth direction. The 
degradation is particularly pronounced in high flying SARs such as the one deployed 
on SEASAT. Azimuthal smearing limits the usefulness of SEASAT data in studying 
the refraction of wave systems across currents. As is discussed in section 4.3, the wave 
system examined in this case study is refracted from range into azimuth. This fact made 
it difficult to track the wave system across the Gulf Stream. 



3.2 Orbit 1339 Images 

The synthetic aperture radar deployed on SEASAT had a lifetime of only about three 
months. There are a limited number of satellite passes which had interesting wind and 
wave conditions along with independent measures of these parameters. Pass 1339 of 
September 28, 1978 was one such pass. 

Beal et al. (1983) and (1986) reported on the processing of the SAR images and the 
oceanographic conditions along the entire pass. The present study focuses on the portion 
of the pass in the vicinity of the Gulf Stream. Figures 3-1 and 3-2 which appeared in 
Beal et al. (1986) show the geographic orientation of pass 1339. 

Digitally processed SAR data from MacDonald Dettwiler and Associates consisting 
of 40 x 40 km images were used for the analysis here. A description of the processing of 
these images is included in Chapter 4. 
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Figure 3-1: SEASAT pass 1339 over the North Atlantic (from Beal et al. 1986). 
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Latitude, N 





Figure 3-2: SEASAT pass 1339 over the Gulf Stream (from Beal et al. 1986). 
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Chapter 4 



Case Study: Gulf Stream 



4.1 Environmental and Meteorological Conditions Dur- 
ing Pass 1339 

SEASAT pass 1339 crossed the Gulf Stream on September 28, 1978 at 1521 GMT. 
Two dominant wave systems existed throughout most of the 900 km pass. The primary 
system consisted of 200 m swell which propagated from southeast to northwest; a sec- 
ondary system of 100 m waves was more variable in direction, but generally traveled 
between a westerly to southwesterly direction. 

Beal et al. (1986) demonstrated that the primary system originated in Tropical 
Depression 12 on September 24. The maximum winds in the vicinity of the depression 
were estimated at 20 m/s and generated waves which traveled towards the west. This 
wave source can be seen in Figure 4-1. 

A smaller storm closer to the satellite pass developed on 26 and 27 September. The 
peak winds were estimated at 15 m/s. This storm can also be seen on Figure 4-1 and is 
the source of the secondary 100 m wave system. 

In the immediate vicinity of the satellite pass, a weak trough dominates the local 
wind field causing a counter clockwise wind movement over the 900 km length of the 
satellite pass. The wind is southerly with a speed of about 5 m/s over the Gulf Stream. 
The wind vectors are shown in Figure 3-1. 
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Figure 4-1: Sources of wave systems observed by SEASAT S AE on pass 1339 (from Beal 
et al. 1986). 
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4.2 SAR Data Analysis 



The available set of SAR data from pass 1339 consists of nineteen 40 x 40 km images 
digitally processed by MacDonald Dettwiler and Associates (Vanvouver, Canada). Beal 
et al. (1983) and (1986) looked at the entire pass to follow the propagation of the 200 
m swell system. In the current study, only the two 40 km images which span the Gulf 
Stream are re-examined in a detailed analysis of the wave-current interactions, as the 
200 m swell system enters the flow of the Gulf Stream at an oblique angle and undergoes 
refraction. 

PVom each of the two images, eleven frames of size 512 x 512 pixels were defined. 
These provided overlapping images of 6.4 x 6.4 km on a side. Each of the twenty-two 
frames thus defined overlaps the next by about 3 km with the exception of the first and 
last in each of the large images. 

The data from each frame is normalized by subtracting the mean value of the frame 
from each pixel value within the frame and then dividing by the mean. A two dimensional 
fast Fourier transform is performed on each of the resulting normalized data frames. The 
complex transform is multiplied by its complex conjugate and oriented with the zero wave 
number value in the center. The resulting transform contains information out to |f 
The image is then tr imm ed to 272 x 272 pixels. An example of a raw Fourier transform 
is shown in Figure 4-2. 

The transforms are spectrally smoothed using a discrete 15 x 15 pixel Gaussian 
function. Following Beal et al. (1986), the smoothed spectrum S, is related to the raw 
spectrum 5„ by: 



f. /» » \ _ Y£?=-7 YltJ-r Su(K+m,k v+ i)e 



(4.1) 



The effect of the smoothing kernel (4.1) on the raw spectrum is also shown in Figure 4- 
2. In this and following figures showing contour plots, adaptive contour intervals were 
used to hi ghligh t energy peaks of interest. The interior of the 256 x 256 pixel smoothed 
region represent wave lengths greater than 50 m. 

The spectra could be further smoothed by along-track averaging of several adjacent 
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(a) 



Unsmoothed Spectrum 




Figure 4-2: (a) Raw FFT of image 20 which is located just after the Gulf Stream, (b) 
Smoothed FFT (same image as (a)). Satellite travel (azimuth direction) is along the 



x-axis. 
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spectra, but this is inappropriate over a region of rapidly changing conditions such as 
the Gulf Stream. 



4.3 Results 

Beal et al. (1986) processed six 6.4 x 6.4 km spectra from each 40 x 40 km SAR 
image. As they followed what appeared to be the dominant 200 m wave system across 
the Gulf Stream, they observed a wave shortening resulting in a wave number increase 
four times that which their model predicted. The wave number vector seemed to turn 
into the current, contrary to the behavior predicted by Snell’s law. Bead et al. concluded 
that it is possible that the system they traced may actually be a secondary system inside 
the Gulf Stream. Figure 4-3 shows the wave length and angle of the system which Beal et 
ad. traced across the Gulf Stream. The wave length behaved approximately as expected 
for a wave which wais propagating into an opposing current at an oblique angle, but the 
direction turned in the wrong sense. 

In this study a closer spacing of wave spectra is used to allow trends to be better 
seen in the rapidly changing environment of the Gulf Stream. It can be seen that merely 
chosing the highest energy spectrad peak does not ensure that the correct system is always 
observed during the wave system’s propagation through the Gulf Stream. Figure 4-4 
show the geographical orientation of the SEASAT SAR images from pass 1339 which 
were processed into the twenty-two spectra. Figures 4-5 through 4-10 display these 
spectra. The energy peak of the 200 m wave system is marked on each figure. It cam be 
seen that azimuthal smearing makes the determination of this peak difficult inside the 
Gulf Stream. The wave lengths and direction of these peaks are plotted on Figure 4- 
11. The shortening and refraction of the wave system toward the North inside the Gulf 
Stream is visible. These results suggest that Beal et al. did apparently trace the wrong 
system through the Gulf Stream. 

It can be seen from Figure 4-11 that the wave length of the wave system across the 
Gulf Stream is shortened from 200 m to about 170 m. At a 40° incidence angle between 
the wave front and the Gulf Stream normal, a 200 m wave should be shortened to 175 m 
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(a) 



Aloog-track distance (km) 




Akjn^track distance fkm) 




Figure 4-3: (a) Wave length and (b) wave angle variation across the Gulf Stream (from 
Beal et aL 1986). 
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Longitude, Degrees West 



Figure 4-4: Geographical orientation of SEAS AT pass 1339 over the Gulf Stream. Num- 
bers at right of track indicate high-resolution SAR image sequence 1 through 22, and 
numbers at left refer to the two larger 40 x 40 km images. 
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Plot 



Plot 3 




Figure 4-5: SAR wave spectra: images 1 through 4. 
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Plot 5 



Plot 7 




Figure 4-6: SAR wave spectra: images 5 through 8. The Gulf Stream edge is located 
between images 7 and 8. Note that the energy peak corresponding to the wave system 
of interest is not evident in image 5. 
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Figure 4-7: SAR wave spectra: images 9 through 12. Note a 6.7 km gap occurs between 
spectra 11 and 12. 
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Plot 13 



Plot 15 





Figure 4-8: SAR wave spectra: images 13 through 16. The North Wall of the Gulf 
Stream is located between images 15 and 16. 
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Plot 17 



Plot 19 




Plot 18 



Plot 20 




Figure 4-9: SAR wave spectra: images 17 through 20. 
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Plot 21 





Figure 4-10: SAH wave spectra: images 21 and 22. 
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Distance along x-axis (m) xlO 4 



Figure 4-11: Wave length and direction of wave front crossing the Gulf Stream. The 
error bars represent the range of possible values along the ridges seen in the spectra in 
the Gulf Stream. The x-axis is oriented along the east-west direction; the waves are 
propagating toward the coast. 
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in the presence of a 2.0 m/s current and to 165 m/s in the presence of a 2.5 m/s current. 
This is in the range of current velocities often observed in the Gulf Stream (Williams 
et al. 1974). Within the Gulf Stream, azimuthal smearing results in the energy center 
appearing as a ridge rather than as a clear peak. The range of possible values along 
this ridge are shown by the error bars within the Gulf Stream on Figure 4-11. The wave 
lengths are more variable than the angles, because the smearing in azimuth degrades the 
wave lengths more than the angles. 

Figures 4-12 and 4-13 show the SAR images reduced in resolution. These allow us 
to visually examine the southern and northern boundaries of the Gulf Stream. The 
southern boundary seen in Figure 4-12 falls at an angle of about 7° or 8° relative to 
the satellite trajectory. Since the satellite pass was 23° from North, this indicates that 
the current was flowing at about 75° from North at the southern boundary of the Gulf 
Stream. The northern boundary is very distinct and can be measured to rise at an angle 
of 7°. This would correspond to a current direction of 60° at the North Wall of the Gulf 
Stream. These measurements indicate that the stream is widening at the point at which 
the satellite passes over. 

The data shown in Figure 4-11 imply that the refraction effect is strongest near the 
northern boundary of the Gulf Stream. This fact would suggest an asymmetric current 
with a velocity peak shifted towards the North Wall. This might also be a reason why 
the boundary is more sharply defined in Figure 4-13 than in Figure 4-12. 

In order to use the SAR data to estimate the current field at the point at which the 
satellite passed over the Gulf Stream, a number of possible currents were examined with 
the refraction model. The width of the Gulf Stream current is estimated to be 33.6 km. 
This is based on measurements from the 40 x 40 km images as well as from the eight 
spectra between the start and end of the Gulf Stream overpass. The spectra have an 
average spacing of 3.36 km, and there is another gap of that size between the 40 x 40 km 
images. The satellite pass direction is nearly normal to the Gulf Stream, so a total of 
ten 3.36 km spacings results in a current approximately 34 km wide. The grid spacing 
used in the numerical model is 51 x 51 grid points each 1.68 km apart (i.e., the model 
grid spacing is half the size of the spectra spacing). The assumption is also made that 
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Figure 4-12: First 40 x 40 km SAR image showing the location of the southern edge of 
the Gulf Stream (arrow). 
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Figure 4-13: Second 40 x 40 km SAR image. The northern boundary (North Wall) of 
the Gulf Stream is clearly visible in the image (arrow). 
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the current direction is constant, so eddies and lateral variations within the Gulf Stream 
are discounted. 

Current directions were proposed to vary from 60° to 75°. The current profiles ex- 
amined vary with a cosine-squared dependence (symmetic jet), a skewed cosine-squared 
dependence with the peak: closer to the northern boundary of the Gulf Stream (asym- 
metric jet), and a top hat shape with tapered edges (top hat). Cross-sections of these 
profiles can be seen in Figure 4-14. Peak current speeds (Umax) within these profiles of 
1.5, 2.0, 2.25 and 2.5 m/s were examined. Each of these current directions, profiles, and 
speeds were input into the numerical model for a 200 m wave with an initial direction 
of 295°. Note that the average angle of the center of the spectra before the Gulf Stream 
was 295°, aind the average wave length was 200 m (11.3 second period). Figures 4-15 
through 4-17 show the orthogonals and rays being refracted across six model currents. 

The wave rays show the propagation of wave energy, so they are compared to the 
spectra processed from the SAR data. Figures 4-18 and 4-19 show ray angles and wave 
lengths from the SAR data (+) superimposed on the model output. Both the model and 
the data are oriented along the x-axis. For each data point, the corresponding model 
values for angle and wave length are found from simple interpolation. From these two 
series of wave lengths and angles, two sets of vectors are defined in the direction of the 
wave ray and with the magnitude of the wave length. The SAR data are compared to the 
model output by finding the phase angle of the complex correlation coefficient between 
the two vector series. This technique is described in Kundu (1975). The vectors are 
defined as: 

u>(t) = u(t ) + iv(t ) (4.2) 

and axe divided into east, u(t), and north, v(t), components. The complex correlation 
coefficient between the model vector series w m (t) and the SAR data vector series a ><f(t) 

, = <ATtW» T ( 4 . 3 ) 

< CJ^(t)w d (t) >3 < a£(t)^m(0 > 2 

where u;’(t) is the complex conjugate. The correlation can also be expressed in terms of 
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east-north components, i.e., 



<*Si + «m> 5 <“3 + »a> 5 

and the phase angle is : 



r + «- 



< t LmVd ~ UdV m > 

<”5, + ''5,>5<“3+<'j>’ 



(4.4) 



, -1 < «mVd - UdVm > 

e avg = tan 

< «m « d + V m V d > 

Note that < ... > indicates the average value of the included series. These complex 
correlations are listed in Tables 4.1, 4.2, and 4.3. All magnitudes of correlation are high 
because the vectors are all pointing close to the same direction. 

To provide another measure of how well the numerical model output fits the SAR 
data, the wave length and angle are examined individually by calculating residuals. The 
residual variance in wave length and angle are calculated from the deviation of the model 
results from the SAR data: 



wave length residual = ^(i u ,/ — m^/) 2 (4.6) 

angle residual = ^(x ang - m ang ) 2 (4.7) 

(4.8) 



where N is the number of data points x w i and x ang are the SAR measured wave lengths 
and angles, and m^i and m ang are the numerical model generated wave lengths and 
angles. 

The complex correlations and residuals are examined together to find the best overall 
fit between the SAR data and the model output in order to determine the current 
characteristics which could explain the observed wave refraction. 

As can be seen from Tables 4.1, 4.2, and 4.3, the best overall correlation between the 
numerical refraction model and the SAR data is found when the current is represented 
by a top hat profile, a velocity of 2.0 m/s, and a direction 75°. In general, the correlation 
improves as the profile is changed from a symmetric jet, to an asymmetric jet, to a top 
hat shape. This indicates that most of the refraction occurs as the wave enters and leaves 
the Gulf Stream, so the current profile within the Gulf Stream at the location of the 
satellite pass is close to uniform, and the largest velocity gradients are near the current 
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(C) 




Figure 4-14: Model current profiles: (a) Symmetric jet. (b) Asymmetric jet . (c) Top 
hat. 
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1.5 m/s symmetric jet current at 65 deg. 




x-axis (m) xlO 4 



xlO 4 2.5 m/s symmetric jet current at 75 deg. 




Figure 4-15: Refraction of wave rays and orthogonals across symmetric jet current pro- 
files. (a) Umax is 1.5 m/s at 65°. (b) U max is 2.5 m/s at 75°. Note that a heavy line 
represents the ray path and a thin line the orthogonal path. 
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1.5 m/s asymmetric jet current at 60 deg 




x 10 4 2.5 m/s asymmetric jet current at 75 deg 




Figure 4-16: Refraction of wave rays and orthogonals across asymmetric jet current 
profiles, (a) U ma T is 1.5 m/s at 60°. (b) C/mo* is 2.5 m/s at 75°. Note that a heavy line 
represents the ray path and a thin line the orthogonal path. 
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1.5 m/s top hat current at 75 deg 




xlO 4 2.0 m/s top hat current at 60 deg 




Figure 4-17: Refraction of wave rays and orthogonals across top hat current profiles, (a) 
Umax is 1.5 m/s at 75°. (b) Umax is 2.0 m/s at 60°. Note that a heavy line represents 
the ray path and a thin line the orthogonal path. 
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(a) 



330 



2.25 m/s symmetric jet current at 75 deg 





Figure 4-18: (a) Symmetric jet profile, (b) Asymmetric jet profile. SAR data (+) 
superimposed on model output; the waves are propagating toward the coast. The x-axis 
corresponds to the east-west direction. 
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Figure 4-19: (a) Top hat profile with 2 m/s current, (b) Top hat profile with 2.25 
m/s current. SAB. data (+) superimposed on model output; the waves are propagating 
toward the coast. The x-axis corresponds to the east-west direction. 
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Umax 


Current 

Direction 


Magnitude of 
Correlation 


Phase 

Angle 


Angle 

Residual deg 2 


Wave Length 
Residual m 2 


1.5 


60° 


.9908 


4.0180° 


89.8 


136.3 


1.5 


65° 


.9916 


3.6014° 


82.1 


125.1 


1.5 


70° 


.9923 


2.9953° 


74.1 


118.2 


1.5 


75° 


.9932 


2.7280° 


67.4 


110.9 


2.0 


60° 


.9922 


3.0911° 


74.6 


127.1 


2.0 


65° 


.9926 


2.3074° 


67.8 


121.8 


2.0 


70° 


.9937 


1.8731° 


60.7 


111.7 


2.0 


75° 


.9952 


1.7376° 


50.5 


113.6 


2.25 


60° 


.9928 


3.0420° 


70.1 


128.8 


2.25 


65° 


.9932 


2.1865° 


63.1 


127.0 


2.25 


70° 


.9943 


1.7043° 


56.7 


114.6 


2.25 


75° 


.9950 


1.1597° 


50.7 


123.8 


2.5 


60° 


.9919 


2.3578° 


73.0 


132.6 


2.5 


65° 


.9928 


1.4241° 


63.0 


133.4 


2.5 


70° 


.9942 


1.3091° 


56.4 


124.8 


2.5 


75° 


.9956 


1.4134° 


48.7 


127.4 



Table 4.1: Correlations of SAE data vs numerical refraction model for symmetric jet 
current profile. 



Umax 


Current 

Direction 


Magnitude of 
Correlation 


Phase 

Angle 


Angle 

Residual deg 2 


Wave Length 
Residual m 2 


1.5 


60° 


.9908 


4.0341° 


89.6 


134.3 


1.5 


65° 


.9918 


3.4259° 


79.7 


121.3 


1.5 


70° 


.9928 


3.0297° 


71.5 


108.9 


1.5 


75° 


.9933 


2.3724° 


64.3 


110.8 


2.0 


60° 


.9927 


3.0574° 


70.4 


123.2 


2.0 


65° 


.9928 


2.0540° 


65.7 


119.2 


2.0 


70° 


.9939 


1.5348° 


57.8 


112.9 


2.0 


75° 


.9958 


1.5629° 


46.8 


106.1 


2.25 


60° 


.9936 


3.1710° 


66.0 


122.1 


2.25 


65° 


.9942 


2.3687° 


58.5 


117.2 


2.25 


70° 


.9953 


1.7400° 


50.2 


110.8 


2.25 


75° 


.9955 


0.6293° 


46.7 


118.3 


2.5 


60° 


.9929 


2.774 7° 


68.2 


126.0 


2.5 


65° 


.9938 


1.4798° 


58.0 


120.0 


2.5 


70° 


.9951 


1.6270° 


51.2 


122.5 


2.5 


75° 


.9962 


2.4065° 


48.4 


122.0 



Table 4.2: Correlations of SAR data vs numerical refraction model for asymmetric jet 
current profile. 
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Umax 


Current 

Direction 


Magnitude of 
Correlation 


Phase 

Angle 


Angle 

Residual deg 2 


Wave Length 
Residual m 2 


1.5 


60° 


.9921 


3.1317° 


75.7 


112.0 


1.5 


65° 


.9940 


2.4323° 


61.7 


95.5 


1.5 


70° 


.9941 


2.1287° 


58.2 


92.8 


1.5 


75° 


.9953 


1.0883° 


49.3 


96.0 


2.0 


60° 


.9930 


1.7350° 


64.3 


108.3 


2.0 


65° 


.9945 


0.9586° 


54.8 


99.7 


2.0 


70° 


.9936 


0.5190° 


48.4 


94.1 


2.0 


75° 


.9963 


0.1357° 


45.1 


98.9 


2.25 


60° 


.9937 


2.0353° 


61.7 


110.2 


2.25 


65° 


.9926 


0.0347° 


66.2 


114.0 


2.25 


70° 


.9965 


1.3584° 


46.4 


106.8 


2.25 


75° 


.9972 


0.5551° 


41.5 


111.2 


2.5 


60° 


.9929 


2.5981° 


70.0 


120.3 


2.5 


65° 


.9943 


2.5297° 


63.4 


119.1 


2.5 


70° 


.9948 


1.8436° 


60.6 


136.0 


2.5 


75° 


.9955 


-.1881° 


54.5 


126.2 



Table 4.3: Correlations of SAR data vs numerical refraction model for top bat current 
profile. 



edges. Tbe correlation improvement from the symmetric jet to the asymmetic jet shape 
indicates that the velocity may be slightly skewed towards the North Wall. 

The current direction which provides the best fit with the SAR, data seems to be 
between 70° and 75°. The southern boundary of the Gulf Stream is oriented at about 
75°, so this value is reasonable. The Umax which provides the best correlation is 2.25 
m/s for the symmetric jet current profile, 2.5 m/s for the asymmetric jet current profile, 
and 2.0 m/s for the top hat profile. 

There is a fair amount of data scatter both before and within the Gulf Stream. This 
makes it diffi cult to assign a specific meaning to a pair of data points which digress from 
the trend followed by the rest. However, both the wave length and angle of propagation 
show two data points within the Gulf Stream which differ significantly from a smooth 
curve. This could be indicative of an eddy formation or lateral variability within the 
Gulf Stream at this point; the length scale is between 6 km and 12 km. 
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Chapter 5 



Summary and Conclusions 



Ocean waves which, interact with currents in deep water undergo refraction, and 
characteristics of the waves like length and height are affected by the interaction. This 
type of refraction can be observed, as swell propagates through major current systems 
like the Gulf Stream or the Antrarctic Circumpolar Current. A model which examines 
wave-current kinematics can be used in conjunction with remote sensing measurements 
of wave propagation to infer the current system which caused the measured refraction. 

The study of the kinematics of wave-current refraction really involves two problems: 
(1) the direction of wave propagation, and (2) the changes in wave length. The wave or- 
thogonals (perpendicular to wave crests) follow paths tangent to the wave number vector 
at a phase speed modified by current interactions. This is the direction of propagation 
seen by an external observer who takes a ’snapshot’ of the wave system. The wave ray 
gives the direction of energy propagation traveling at the group velocity and is given by 
the vector sum of the intrinsic group velocity and the current velocity. 

A numerical model is developed which predicts the direction of propagation and wave 
length, as a wave ray or orthogonal is refracted by a current. This model is validated by 
comparing its results to analytical solutions for cases of simple refraction problems. 

SEASAT SAR data from a pass over the Gulf Stream are analyzed to allow wave 
energy propagation to be traced across the Gulf Stream. The digitally processed SAR 
images are sectioned into overlapping 512 x 512 pixel frames. These frames are normal- 
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ized, and a Fast Fouiier Transform is performed on each one. The magnitude of the 
transform is found, and it is oriented with the zero wave number at the center. The 
resulting raw spectra are smoothed and trimmed to obtain energy spectra which show 
waves of length greater than 50 m. These spectra contain energy peaks which corre- 
spond to a 200 m wave system which traveled from southeast to northwest across the 
Gulf Stream. The peaks are distinctly visible before and after the Gulf Stream, but 
within the Gulf Stream the waves are refracted into azimuth, and the effect of azimuth 
or Doppler smearing makes it much more difficult to trace the wave energy propagation. 

The general position and width of the Gulf Stream below the satellite pass is clearly 
marked in the SAR images. Three possible current profiles are proposed to model the 
Gulf Stream: a symmetric jet, an assymetric jet, and a top hat. The current direction is 
modeled by directions of flow which vary from 60° to 75°, and the current velocity, Umax, 
is varied from 1.5 to 2.5 m/s. These current parameters are input into the numerical 
model as a current field. Using the model, the kinematics of a 200 m wave passing 
through the Gulf Stream are examined. The model output is compared to the SAR data 
by colocating the model with the data by representing the wave length and direction of 
propagation as sets of vectors. From this the magnitude and phase angle of the complex 
correlation between the SAR data and the model are calculated. The wave length and 
angle residuals are also calculated to provide an additional measure for comparison. 

The best correlations are found for Gulf Stream velocities between 2.0 and 2.5 m/s 
at 70° to 75°. The best overall correlation between the model and SAR data is obtained 
when the Gulf Stream is modeled by a top hat profile, a velocity of 2.0 m/s, and a 
direction of 75°. 75° agrees with the direction visually observed from the SAR images, 
and the direction and speeds are close to the Coast Guard estimates for the Gulf Stream 
on this day. The current profiles used in the model do not take into account a possible 
widening of the Gulf Stream at the position of the satellite overpass. There is a great 
deal of scatter in the SAR data both before and in the Gulf Stream, so it is difficult to 
correlate every point with specific current behavior, but the increase in wave length and 
change in angle in the center of the Gulf Stream seem to indicate that there may be a 
non-uniform feature such as an eddy or lateral variation forming within the current. 
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SAR data can be used to trace wave systems which pass through a current like the 
Gulf Stream. The tomographic employment of SAR data in conjunction with a numerical 
model for wave-current refraction can be used to infer the current which caused the wave 
refraction. In order for the approach described above to be practically useful, a uniform 
and stable long wave system must exist. This tomographic application requires a wave 
system which stays in the range direction with regard to the satellite, or which is long 
enough to still be measured in azimuth. The large scale current behavior is quite evident 
in the SAR data, but small scale eddies or lateral variations are harder to measure, 
because they may have smaller length scales than the scales of the frames processed 
here. 

Lower flying SARs which overcome some of the inherent problems of azimuth smear- 
ing sire superior to the SAR which flew in SEASAT for this type of tomographic ap- 
plication. Two Shuttle missions deployed SARs, the SLR-A and SIR-B flights; they 
provided relatively short sets of SAR data from lower altitudes than SEASAT. Even 
SARs deployed in high flying aircraft can provide wide geographic coverage of current 
systems. 

SAR and numerical models of wave-current refraction offer promising new tools to 
the oceanographer, as he seeks to better measure and understand the ocean environment. 
When they are used together, the structure of current systems can be measured from 
space. This type of tomography is limited by azimuth smearing, as a wave system is 
refracted out of the range and into the azimuth direction. Lower flying SARs which will 
overcome azimuth smearing and also provide better small scale resolution are needed to 
continue this area of oceanographic research. 
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